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Abstract. 

We analyse the effect of the compression of the dark matter due to the infall of baryons to the galactic center on the 
gamma-ray flux. In addition, we also consider the effect of non-universal supersymmetric soft terms. This analysis shows 
that neutralino dark matter annihilation can give rise to signals largely reachable by future experiments like GLAST. This is 
a remarkable result if we realise that direct detection experiments will only be able to cover a small region of the parameter 
space. Actually, in this SUGRA framework we have also been able to fit present excess from EGRET and CANGAROO using 
different non-universal scenarios, and even fit the data from both experiments with only one scenario. We have also studied 
the recent HESS data implying a neutralino heavier than 12 TeV. Because of such a heavy neutralino, it is not natural to find 
solutions in the SUGRA framework. Nevertheless we have carried out a quite model-independent analysis, and found the 
conditions required on the particle physics side to fit the HESS data thanks to dark matter annihilation. 
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INTRODUCTION 

It is now well established that luminous matter makes 
up only a small fraction of the mass observed in Uni- 
verse. A weakly interacting massive particle (WIMP) is 
one of the leading candidates for the "dark" component 
of the Universe. One of the most promising methods for 
the indirect detection of WIMPs consists of detecting the 
gamma rays produced by their annihilations in the galac- 
tic halo. Concerning the nature of WIMPs, the best mo- 
tivated candidate is the lightest neutralino, a particle pre- 
dicted by the supersymmetric (SUSY) extension of the 
standard model (lj. We implement in our analysis |01 
the lower bounds on the masses of SUSY particles and 
Higgs boson, as well as the experimental bounds on the 
branching ratio of the b — > sy process and on a^ USY , for 
which the more stringent constraint from e + e~ disfavors 
important regions of the SUGRA parameter space (see 
e.g. Ref. 0,01). In addition, we have also taken into ac- 
count the last data concerning the B s — > branching 
ratio. 

THE GAMMA-RAY FLUX 

For the continuum of gamma rays, the observed differen- 
tial flux at the Earth coming from a direction forming an 



angle \jf with respect to the galactic center is 
dN' 1 r 
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where the discrete sum is over all dark matter annihi- 
lation channels, dN'y/dEy is the differential gamma-ray 
yield, (<7,v) is the annihilation cross section averaged 
over its velocity distribution, m x is the mass of the dark 
matter particle, and p is the dark matter density. Ac- 
tually, when comparing to experimental data one must 
consider the integral of 4> r over the spherical region of 
solid angle Ai2 given by the angular acceptance of the 
detector which is pointing towards the galactic center. 
For example, for EGRET A£2 is about 10~ 3 sr whereas 
for GLAST, CANGAROO, and HESS it is 1(T 5 sr. Note 
that neutralinos move at galactic velocity and therefore 
their annihilation occurs at rest. 



THE ADIABATIC COMPRESSION 
MODEL 

Recently, SUSY dark matter candidates have been stud- 
ied in the context of realistic halo models including bary- 
onic matter |5]. Indeed, since the total mass of the inner 
galaxy is dominated by baryons, the dark matter distri- 
bution is likely to have been influenced by the baryonic 
potential. In particular, its density is increased, and as a 
consequence typical halo profiles such as Navarro, Frenk 



and White (NFW) |0] and Moore et al. 0] have a more 
singular behaviour near the galactic center. The conclu- 
sion of the work in Ref. |5] is that the gamma-ray flux 
produced by the annihilation of neutralinos in the galac- 
tic center is increased significantly, and is within the sen- 
sitivity of incoming experiments, when density profiles 
with baryonic compression are taken into account. In- 
deed, highly cusped profiles are deduced from N-body 
simulations. In particular, NFW |6] obtained a profile 
with a behaviour p(r) °< r _1 at small distances. A more 
singular behaviour, p (r) °< r~ 15 , was obtained by Moore 
et al. 13]. However, these predictions are valid only for 
halos without baryons. One can improve simulations in 
a more realistic way by taking into account the effect of 
the normal gas (baryons). This loses its energy through 
radiative processes falling to the central region of form- 
ing galaxy. As a consequence of this redistribution of 
mass, the resulting gravitational potential is deeper, and 
the dark matter must move closer to the center increasing 
its density. 

This increase in the dark matter density is often treated 
using adiabatic invariants. The present form of the adi- 
abatic compression model was numerically and analyti- 
cally studied by Blumental et al. 1 8 ] . This model assumes 
spherical symmetry, circular orbit for the particles, and 
conservation of the angular momentum M(r)r = const., 
where M(r) is the total mass enclosed within radius r. 
The mass distributions in the initial and final config- 
urations are therefore related by M,(r,)r/ = [M&(r/) + 
MoAi(rf)]rf, where M,(r), M b (r) and M DM (r) are the 
mass profile of the galactic halo before the cooling of 
the baryons (obtained through N-body simulations), the 
baryonic composition of the Milky Way observed now, 
and the to be determined dark matter component of the 
halo today, respectively. This approximation was tested 
in numerical simulations |9]. Nevertheless, a more pre- 
cise approximation can be obtained including the possi- 
bility of elongated orbits |5]. The models and constraints 
that we used in this work for the Milky Way can be found 
in Table I of Ref. |5J. As one can see in |2J], at small r the 
dark matter density profile following the adiabatic cool- 
ing of the baryonic fraction is a steep power law p oc r^ c 
with y c w 1.45(1.65) for a NFW c (Moore c ) compressed 
model. 

SUPERSYMMETRIC MODELS 

As discussed in detail in Ref. 0] in the context of indirect 
detection, <7,- can be increased in different ways when the 
structure of mSUGRA for the soft terms is abandoned. 
In particular, it is possible to enhance the annihilation 
channels involving exchange of the CP-odd Higgs, A, 
by reducing the Higgs mass. In addition, it is also pos- 
sible to increase the Higgsino components of the lightest 



neutralino. Thus annihilation channels through Higgs ex- 
change become more important than in mSUGRA. This 
is also the case for Z T , xf, and j^-exchange channels. 
As a consequence, the gamma-ray flux will be increased. 

In particular, the most important effects are produced 
by the non-universality of Higgs and gaugino masses. 
These can be parameterised, at Mqut, as follows 

m^=m 2 (l+5 1 ), m 2 Hu =m 2 {\ + & 2 ) , (2) 

and 

M,=M, M 2 =M(l + 8£) , M 3 =M(l + 8i), (3) 
We concentrated in [2] our analysis on the following 



representative cases: 
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Case a) corresponds to mSUGRA with universal soft 
terms, cases b), c) and d) correspond to non-universal 
Higgs masses, and finally cases e) and/) to non-universal 
gaugino masses. The cases b), c), d), and e) were dis- 
cussed in Ref. J4J], and are known to produce gamma-ray 
fluxes larger than in mSUGRA, whereas case /) will be 
of interest when discussing heavy WIMP signals predic- 
tions in the perspective of atmospheric Cherenkov tele- 
scopes like e.g. CANGAROO. 

CONFRONTING EXPERIMENTS 
EGRET 

The EGRET telescope on board of the Compton 
Gamma-Ray Observatory has carried out the first all- 
sky survey in high-energy gamma-rays (« 30 MeV - 
30 GeV) over a period of 5 years, from April 1991 un- 
til September 1996. As a result of this survey, it has de- 
tected a signal H 1 Oil above about 1 GeV, with a value for 
the flux of about 10~ 8 cm~ 2 s~\ that apparently can- 
not be explained with the usual gamma-ray background. 
The source, possibly diffuse rather than pointlike, is lo- 
cated within the 1 .5" (AO. ~ 10~ 3 sr) of the galactic cen- 
ter. Due to the lack of precision data in the high energy 
bins, it seems impossible however to distinguish any an- 
nihilation channel leading to this photon excess. The re- 
sults can be seen in Fig. [2 where case c) of Eq. (0} have 
been studied for a scan on m and M from to 2 TeV 
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FIGURE 1. Gamma-ray spectra <fry(Ey) from the galactic 
center as functions of the photon energy for the SUGRA case 
c) discussed in Eq. J4j for tan/3 = 35, A = and fl > 0, 
compared with data from the EGRET experiment. NFW and 
Moore et al. profiles with adiabatic compression are used with 
Afl ~ 1CP 3 sr. All points shown after a scan on m and M (0 
- 2000 GeV) fulfil the accelerator constraints discussed in the 
text, and WMAP bounds. 



and tan/3 = 35. Let us remark that it is possible to differ- 
entiate each point of the parameter space (to, M) by its 
gamma-ray spectrum. The higher fluxes for instance cor- 
responds to the closing of the A-pole, whereas the lower 
flux spectrum is obtained through the opening of the A- 
pole (see e.g. point B of Fig. 3 in |yl). 



CANGAROO-GLAST 

Recently, the CANGAROO-II atmospheric 
Cherenkov telescope has made a significant detec- 
tion of gamma rays from the Galactic center region 
lITlll . In particular, the collaboration has published the 
spectrum obtained in six energy bins, from 200 GeV to 
3 TeV. Observations taken during 2001 and 2002 have 
detected a statistically significant excess at energies 
greater than ~ 250 GeV, with an integrated flux of 
~2x 10~ 10 photons cm~ 2 s _1 . These measurements 
indicate a very soft spectrum °< /j- 4 - 6 * - 5 . 

It is interesting to see whether it is possible to obtain 
such a candidate in SUGRA scenarios imposing the ac- 
celerator and WMAP constraints, and withing the frame- 
work of adiabatically compressed halos. The baryonic 
cooling effect on the fluxes gives us the order of magni- 
tude needed to fit with both data with a 1 TeV neutralino 
in the non-universal case e) with M3 = 0.5M. It is worth 
noticing that the CANGAROO-II collaboration in (TT1 
pointed out already that the EGRET and CANGAROO- 
II data can be relatively smoothly connected with a cutoff 
energy of 1-3 TeV. Typical points of the parameter space 
fullfilling all experimental constraints and fitting both set 
of data lie between (to = 800 GeV, M = 800 GeV) and 
(to = 3 TeV, M = 3 TeV). 
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FIGURE 2. Gamma-ray spectra from the galactic center as 
functions of the photon energy for the non-universal case 
e) with M3 = 0.5M, compared with data for EGRET and 
CANGAROO-II experiments and the expected GLAST sen- 
sitivity. Here only the average profile defined in Sect. ?? with 
adiabatic compression, NFW^., is used. 



It is also interesting to see the complementarity of 
GLAST with EGRET and CANGAROO. GLAST will 
perform an all-sky survey detection of fluxes with energy 
from 1 GeV to 300 GeV, exactly filling the actual lack 
of experimental data in this energy range (see Fig. |5Jl, 
and checking the CANGAROO results. Indeed, we have 
calculated that the integrated gamma ray flux for such 
a signal will be around 5 x 10~ n cm~ 2 s _1 . We have 
shown this sensitivity curve in Fig. [2] for AQ. = 10~ 5 , 
which is the typical detector acceptance, following the 
prescriptions given in 11311 . We clearly see that GLAST 
will finish to cover the entire spectrum. 



HESS 

The HESS Cherenkov telescope experiment has re- 
cently published new data on gamma rays, detecting a 
signal from the Galactic Center 11211 . The measured flux 
and spectrum differ substantially from previous results, 
in particular those reported by the CANGAROO collabo- 
ration, exhibiting a much harder power-law energy spec- 
trum with spectral index of about —2.2 and extended up 
to 9 TeV. The authors of 1 12] already pointed out that 
if we assume that the observed gamma rays represent 
a continuum annihilation spectrum, we expect m x > 12 
TeV. Actually such a heavy neutralino-LSP is not natu- 
ral in the framework of a consistent supergravity model 
when we impose the renormalisation group equations 
and radiative electroweak symmetry breaking. 

Although in |2] we performed some scans in all non 
universality directions using numerical dichotomy meth- 
ods, no point in the parameter space in any non-universal 
case studied was able to give a several 10 TeV neutralino 
satisfying WMAP constraint but this can be sensitive to 
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FIGURE 3. Dark matter annihilation versus HESS data for 
NFW compressed halo models. 



the RGE and relic density calculation codes. 

Nevertheless, without RGE and taking soft parame- 
ters at the electroweak scale, the constraints are easier to 
evade. On top of that, in a very effective approach using 
completely free parameters and couplings in cross sec- 
tions, neutralinos with m z ~ 10 TeV and £l x h 2 ~ WMAP 
may certainly be fine tuned. We did not adopt such ap- 
proaches since the MSSM is motivated by high energy 
and theoretical considerations. We analyzed in a quite 
model-independant way the conditions required on the 
particle physics field to fit with the HESS data thanks to 
dark matter annihilation. The results are shown in Fig. [3] 
We point the fact that we do not only compare the spec- 
trum shape of the signal with possible dark matter an- 
nihilation explanation. Indeed, it should be noticed that 
our compressed halo profiles give rise to absolute gamma 
fluxes within the HESS data order of magnitude with 
(ffv) values in agreement with the WMAP requirement. 
It is worth noticing here than a similiar and more com- 
plete analysis was done recently by S. Profumo in 1 14]. 

CONCLUSION 

We have analysed the effect of the compression of the 
dark matter due to the infall of baryons to the galactic 
center on the gamma-ray flux. In addition, we have also 
consider the effect of non-universal soft terms, that arises 
naturally in string motivated framework 1 1 5 ] . This anal- 
ysis shows that neutralino dark matter annihilation can 
give rise to signals largely reachable by future experi- 
ments like GLAST. This is a remarkable result if we re- 
alise that direct detection experiments will only be able 
to cover a small region of the parameter space. Actually, 
in this SUGRA framework we have also been able to 
fit present excess from EGRET and CANGAROO using 
different non-universal scenarios, and even fit the data 
from both experiments with only one scenario. We have 
also studied the recent HESS data implying a neutralino 



heavier than 12 TeV. We have carried out a quite model- 
independent analysis, and found the conditions required 
on the particle physics side to fit the HESS data thanks 
to dark matter annihilation. 

In any case, we must keep in mind that the current data 
obtained by the different gamma-rays observations from 
the Galactic Center region by these three experiments do 
not allow us to conclude about a dark matter annihilation 
origin rather than other, less exotic astrophysics sources. 
Fortunately, this situation may change with the improve- 
ment of angular resolution and energy sensitivity of fu- 
ture detectors like GLAST, which will be a complemen- 
tary source of gamma-ray data. 

ACKNOWLEDGMENTS 

Y.M. want to thank the organisation commity for having 
made him discover the true mystery of Kimchis. 

REFERENCES 

1. For a recent review, see C. Munoz, Int. J. Mod. Phys. A19 
(2004) 3093 | arXiv:hep-ph/0309346 1. 

2. Y. Mambrini, C. Munoz, E. Nezri and F. Prada, 
arXiv:hep-ph/0506204 

3. G. Bertone, E. Nezri, J. Orloff and J. Silk, Phys. Rev. D70, 
063503 (2004) |arXiv:astro-ph/0403322 |. 

4. Y. Mambrini and C. Munoz, arXiv:hep-ph/0407158 
published in Astropart. Phys.; 

Y. Mambrini and C. Munoz, JCAP 10 (2004) 003 
|arXiv:hep-ph/0407352|. 

5. F. Prada, A. Klypin, J. Flix, M. Martinez and E. Simonneau, 
arXiv:astro-ph/0401512 

6. J. F. Navarro, C. S. Frenk and S. D. M. White, Astrophys. 
J. 462 (1996) 563 |arXiv:astro-ph/9508025 |. 

7. B. Moore, T. Quinn, F. Governato, J. Stadel and 

G. Lake. Mon, Not. Roy. Astron. Soc. 310 (1999) 1147 
| arXiv : astro-ph/9903 1 64 1 . 

8. G. R. Blumenthal, S. M. Faber, R. Flores and J. R. Primack, 
Astrophys. J. 301 (1986) 27. 

9. O. Y. Gnedin, A. V. Kravtsov, A. A. Klypin and D. Nagai, 
Astrophys. J. 616 (2004) 16 |arXiv:astro-ph/0406247 |. 

10. EGRET Collaboration, S. D. Hunger et al., Astrophys. J. 
481 (1997) 205; H. A. Mayer-Hasselwander et al., Astron. 
& Astrophys. 335 (1998) 161. 

11. CANGAROO-II Collaboration, K. Tsuchiya et al, 
Astrophys. J. 606 (2004) LI 15 | arXiv:astro-ph/0403592 1. 

12. HESS Collaboration, F. Aharonian et al., Astron. & 
Astrophys. L13 (2004) 425 | arXiv:astro-ph/0408145 |. 

13. N. Fornengo, L. Pieri and S. Scopel, Phys. Rev. D70 
(2004) 103529 |arXiv:hep-ph/0407342 |. 

14. S. Profumo, arXiv:astro-ph/05086^y 

15. P. Binetruy, Y. Mambrini and E. Nezri, Astropart. Phys. 
22 (2004) 1 | arXiv:hep-ph/0312155 |; 

G. Bertone, P. Binetruy, Y. Mambrini and E. Nezri, 
Astropart. Phys. 24 (2005) 44 | arXiv:hep-ph/0406083 |. 
P. Binetruy, A. Birkedal-Hansen, Y. Mambrini and 
B. D. Nelson, arXiv:hep-ph/0308047 



